MULTITERMINAL COUNTING STATISTICS 



DMITRI A. BAGRETS, YULI V. NAZAROV 
Department of Applied Physics and 

Delft Institute of Microelectronics and Suhmicrontechnology, 
Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, 
The Netherlands 



Abstract. The review is given of the calculational schemes that aUows for easy 
evaluation of full current statistics (FCS) in multi-terminal mesoscopic systems. 
First, the scattering approach by Levitov et.al to FCS is outlined. Then the multi- 
terminal FCS of the non-interacting electrons is considered. We show, that this 
theory appears to be a circuit theory of 2 x 2 matrices associated with Keldysh 
Green functions. Further on the FCS in the opposite situation of mesoscopic sys- 
tems placed in a strong Coulomb blockade limit is discussed. We prove that the 
theory of FCS in this case turns out to be an elegant extension of the master 
equation approach. We illustrate both methods by applying them to the various 
three-terminal circuits. We study the FCS of electron transfer in the three-terminal 
chaotic quantum dot and compare it with the statistics of charge transfer in the 
Coulomb blockade island with three leads attached. We demonstrate that Coulomb 
interaction suppresses the relative probabilities of big current fluctuations. Finally, 
for the generic case of single resonance level the equivalence of scattering and mas- 
ter equation approach to FCS is established. 



1. Introduction 

The field of quantum noise in mesoscopic systems has been exploded during 
the last decade, some achievements being summarized in a recent review 
article. [1] While in classical systems shot noise is just a straightforward 
manifestation of discreteness of the electron charge, it can be used in quan- 
tum system as unique tool to reveal the information about the electron 
correlations and entanglement of various kinds. Measurement of fractional 
charge in Quantum Hall regime [2], noise measurements in atomic-size 
junctions [3], chaotic quantum dots[4] and superconductors [5, 6] are mile- 
stones of the field of quantum noise. Starting from the pioneering work of 
Biittiker[7], a special attention has been also paid to shot noise and noise 
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correlations in the "multi-terminal" circuits. The correlations of currents 
flowing to/from different terminals can, for instance, reveal Fermi statistics 
of electrons. Such cross-correlations have been recently seen experimentally 
in Ref. [8, 9, 10]. 

A very important step in the field of quantum noise has been made by 
Levitov et. al. in [11, 12] where an elegant theory of full counting statistics 
(FCS) has been presented. This theory provides an efficient way to investi- 
gate the current correlations in the mesoscopic systems. The FCS method 
concentrates on the evaluation of the probability distribution function of 
the numbers of electrons transferred to the given terminals during the given 
period of time. It yields not only the noise, but all higher momenta of the 
charge transfer. The probability distribution also contains the fundamental 
information about big current fluctuations in the system. 

Original FCS method [11, 12] was formulated for the scattering ap- 
proach to mesoscopic transport. This made possible the study the statis- 
tics of the transport through the disordered metallic conductor [13] and 
the two-terminal chaotic cavity [14]. Muzykantskii and Khmelnitskii have 
generalized the original approach by Levitov et. al. to the case of the nor- 
mal metal/superconducting contacts. The very recent application of the 
scattering approach is the counting statistics of the charge pumping in the 
open quantum dots [15, 16, 17]. 

However, scattering approach becomes unpractical in case of realistic 
layouts where the scattering matrix is multi-channel, random and cum- 
bersome. For practical calculations one evaluates the FCS with the semi- 
classical Keldysh Green function method [18] or with its reduction called 
the circuit theory of mesoscopic transport [19]. The Keldysh method to 
FCS was first proposed by one of the authors. It has been applied to treat 
the effects of the weak localization corrections onto the FCS in the dis- 
ordered metallic wires and later on to study the FCS in superconducting 
heterostructures [20]. 

The above researches address the FCS of the non-interacting electrons 
in case of common two-terminal geometry. Although the cross-correlations 
for several multi-terminal layouts have been understood [1], the evaluation 
of FCS encountered difficulties. For instance, an attempt to build up FCS 
with "minimal correlation approach" [21] has lead to contradictions [14]. 
Very few is known about FCS of interacting electrons [22, 23]. Since the 
interaction may bring correlations and entanglement of electron states the 
study of FCS of interacting electrons, particularly for the case of multi- 
terminal geometry, is both challenging and interesting. 

In the present work we review the calculational scheme that allows for 
easy evaluation of FCS in multi-terminal mesoscopic systems. The paper is 
organized as follows. In the section II we outline the scattering approach by 
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Figure 1. The "multi-terminal" mesoscopic conductor. 



Levitov et.al to FCS. Section III is devoted to the multi-terminal FCS of 
the non-interacting electrons. We show, that this theory appears to be a cir- 
cuit theory of 2 x 2 matrices associated with Keldysh Green functions [24] . 
In the section IV we concentrate on FCS in the opposite situation of meso- 
scopic systems placed in a strong Coulomb blockade limit. We prove that 
the FCS in this case turns out to be an elegant extension of the master 
equation approach. We illustrate both methods by applying them to the 
various three-terminal circuits. We study the FCS of electron transfer in the 
three-terminal chaotic quantum dot and compare it with the statistics of 
charge transfer in the Coulomb blockade island with three leads attached. 
We demonstrate that Coulomb interaction suppresses the relative probabil- 
ities of big current fluctuations. In section V we establish the equivalence 
of scattering and master equation approach to FCS, by considering the 
statistics of charge transfer through the single resonance level. Finally, we 
summarize the results in the section VI. 

2. Scattering approach to FCS 

In this section we review the current statistics of non-interacting electrons, 
using the ideas of general scattering approach to mesoscopic transport. The 
problem of FCS in this framework has been solved about a decade ago by 
Levitov et. al [11, 12]. Their work, concentrated on the two-terminal geome- 
try [12], became widely cited. But, unfortunately, their preceding paper[ll], 
devoted to multi-terminal systems as well, is only partially known. We will 
also demonstrate its relation to the current-current correlations in multi- 
terminal conductors, first investigated by M.Biittiker [7]. 
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The most general mesoscopic system, eligible to scattering approach, is 
shown in Fig. 1. It is a phase coherent mesoscopic conductor, connected 
to macroscopic reservoirs (leads, terminals) by means of n external con- 
tacts (n > 2). We concern the situation when the system is placed in the 
non-equilibrium condition, thus each terminal a is characterized by a dis- 
tribution fTinction fa{£) of electrons over the energy. When the terminal is 
in local equilibrium at electro-chemical potential Va and a temperature Tq, , 
then fa{£) reduces to Fermi distribution /^(e) = {l+ex.p[{e—eVa)/kTa]}~^. 

The conductor is a region of disordered or chaotic scattering. The scat- 
tering within conductor is assumed to be elastic. This is true at sufficiently 
low energies, such that the energy dependent length of inelastic scattering 
exceeds the sample size. As to the terminals, it is conveniently to assume, 
that far from the contacts the longitudinal (along the lead) and transverse 
(across the lead) motion of electrons arc separable. Then at fixed energy e 
the transverse motion in the lead a is quantized and described by the integer 
index Ua = ^ ■ ■ ■ Maie), with Mais) being the total number of propagating 
modes at given energy. Each n„ corresponds to one incoming to the contact 
and one outgoing from the contact a sohition of the Shrodinger equation, 
which is usually referred scattering channel. 

The idea behind the scattering approach is that the electron trans- 
port can be described as the one-particle scattering from the incoming 
channels to the outgoing channel in the leads. Very generally it can be 
described by a unitary scattering matrix S{e). At given energy its dimen- 
sion is equal to M , M = being the total number of transport 
channels. The energy dependence of the scattering matrix is set by the 
"Thouless energy" E^h = h/rd, Td being the traversal time through the 
system. At sufficiently low temperatures, voltages and frequencies, such 
that m.SLx{eV,kT,hLo} <C E^h, one can conveniently disregard the energy 
dependence and concentrate on the scattering matrix S(Ep), taken at Fermi 
surface. Then the scattering approach gives the adequate description of the 
transport, provided the effects of Coulomb blockade can be neglected. This 
is justified, when the total conductance G of the sample is much larger that 
the conductance quantum Gq = e'^/2'Kh. 

One starts by introducing the scattering matrix S of electrons at a 
narrow energy strip ~ inax{eV, kT) <C -Exh in the vicinity of Fermi energy 
Ep- The ^-matrix relates the amplitudes ip^^p of incoming electrons to the 
amplitudes tpii^a of the outgoing ones (See Fig. 1): 




(1) 



m,/3 



Here Latin indexes (n, m) refer to the particular quantum scattering chan- 
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nel in given leads {a, (See e.g. [1, 7, 30] for a more detailed introduction.) 

It is worth to mention that at this stage ipr^l correspond to the amplitude 
of the wave function in the exterior area of the contacts. Thus defined 
S-matrix incorporates the information about scattering in the whole sys- 
tem. This includes the scattering from the disordered/chaotic regions of the 
conductor as well as multiple reflections from external contacts. 

Remarkably, the FCS can be expressed in terms of the scattering matrix 
in quite a general way. It has been shown by Levitov et. al. [11]. The FCS 
approach deals with the probability distribution P{{Ni}) for Ni electrons to 
be transferred through the terminal i during time interval Iq. The detection 
time to is assumed to be much greater than e/I, I being a typical current 
through the conductor. This ensures that in average iVj 3> 1. The prob- 
ability distribution P{{Ni}) can be conveniently expressed via generating 
function ("action") S{{xi}) by means of Fourier transform 

P({iV,}) = r n ^e-'^toiW-E.iv.x.. (2) 

Here Xi ^'^^ the parameters of this generating function ("counting fields"), 
each of them being associated with a given terminal i. The generating func- 
tion S{{xi}) contains the whole information about the irreducible moments 
of charge transfer through the system, as well as the information about big 
current fluctuations. Indeed, it follows from Exp. (2) that (higher-order) 
derivatives of S with respect to Xii evaluated at Xi = 0; gi'^e (higher-order) 
irreducible moments of P{{Ni\). First derivatives yield average currents to 
terminals, second derivatives correspond to the noises and noise correla- 
tions. 

The result of Levitov et. al. for the generating function S{{xi}) reads 
as follows. One considers the energy-dependent determinant 

SE{{Xi}) = lnDet(l -fE + IeS^s) (3) 

where the matrix {fE)rnn,ai3 = fa{E) is diagonal in the channel indexes 
{m,n) and the matrix S is defined as Smn,ai3 = Smn,a(}- It has 

been proved that this expression represents a characteristic function of the 
probability distribution of transmitted charge for electrons with energies in 
an infinitesimally narrow energy strip in the vicinity of E. After that the 
complete generating function S{{xi}) is just a sum over energies 

S{{Xi})=t^ I SEiiXr})^ (4) 

This expression refiects the fact that electrons at different energies are 
transferred independently, without interference. Therefore they yield addi- 
tive independent contributions to the generating function. 
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The general results for the FCS, contained in Eqs. (3) and (4), provides 
the elegant way to derive the pioneering results of M. Biittiker [7], concern- 
ing the current-current correlations in the multiterminal conductors. One 
considers the correlation function 

Papit -t') = l (^AUt)Aifs{t') + Aip{t')AUt)) (5) 

of currents fluctuations in contacts a and /3. The fluctuation AIa{t) is 
defined as AIa{t) = Ia{t)~{Ia), with Ia{t) being the current operator in the 
lead a and (Iq,) being its mean steady value under given non-equilibrium 
conditions. The Fourier transform PQ^(a;) of current-current correlations 
functions is sometimes referred to as noise power. Throughout this article 
we concentrate on zero-frequency (shot-noise) limit of current correlations 
Pap = Pafii^ = 0)- At low frequencies both (7^) and Pap are readily 
expressed via the first (second) moments of number of transferred electrons 
through the corresponding leads 

2 

{la) = f (iV„), Pap = ^{ANaANp) (6) 

to To 

Therefore the second derivative d^^^^'S{{Xi}) of the action yields the re- 
sult for the shot-noise correlation function P^p. It can be reduced to the 
following form 

Pa(3 = ^[ dEj^T^ {^^Al} - Um (7) 

Here the matrix A^,, is defined via 5-matrix as 

A^^{E) = Sa^Sa^, - Si^{E)Sa^,{E) (8) 

This is precisely the result obtained by M. Biittiker prior the development 

of FCS approach. 

Up to the moment the scattering matrix S was not specified. In order 
to evaluate the FCS, one needs to construct such a matrix, so that to take 
into account the scattering properties of the concrete mesoscopic system at 
hand. To accomplish this task, one may proceed as follows. Very generally, 
one can separate a mesoscopic layout into primitive elements: nodes and 
connectors. The nodes are similar to the terminals, in the sense, that each 
of them can be characterized by some non-equilibrium isotropic distribu- 
tion function of electrons over the energies. The illustrative example of the 
node is chaotic quantum dot. The scattering within the node is assumed 
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Figure 2. The building blocks of a scattering matrix of the multi-terminal chaotic 
quantum dot. The matrix Sa describes the scattering due to the contact a. The matrix 
S''''-' describes the chaotic scattering within the dot. 



to be random (chaotic). The connectors represent themselves either inter- 
nal or external contacts in the systems, scattering properties of which are 
supposed to be known. The actual implementation of connectors can be 
quite different: quantum point contacts, tunnel junctions, diffusive wires, 
etc. The purpose of the above separation is to ascribe to each node and to 
connector the particular unitary scattering matrix with known properties. 
Then the scattering matrix S of the entire system can be unambiguously 
expressed via these primitive scattering matrices. 

We outline this scheme for the simple system, consisting of the sin- 
gle node only. A good example of such system is a multi-terminal chaotic 
quantum dot, shown in Fig. 2. It is coupled to n external reservoirs via 
n contacts. Each contact a is described by a unitary scattering matrix 
Sa. (Its choice depends on the experimental realization of this given con- 
tact.) Let 4'm}a be the coefficients of electron wave function, taken at the 

interior boundary of the contact. The amplitudes (t)m}a and (t)m}a cor- 
respond to the outgoing (incoming) waves from (to) the the contact a 
respectively, as shown in Fig. 2. The scattering matrix Sa relates the 
vector c™ = {'4'Q^\4'^a ^)^ of the incident electron waves to the vector 
gOut _ ^^(^ ) ^ 4>'a^)'^ of the outgoing waves: 

c"^' = Sad^=(^ '^d- (9) 

Here the square blocks (r^) of the size {Ma x Ma) describe the elec- 
tron reflection back to the reservoirs or back to the dot, respectively. The 



8 DMITRI A. BAGRETS, YULI V. NAZAROV 

off-diagonal blocks ta (i^) of the same dimension are responsible for the 
transmission through the contact a. 

Similarly, the chaotic scattering of electrons inside the dots is described 
by the scattering matrix S^^\ By analogy to Eq. (1), it relates within the 

dot amplitudes ^^^^ to ^n, a, i-e. 

The dimension (M x M) of S^^^ is equal to that of the matrix S. 

It is possible now to express the scattering matrix S of the entire sample 
in terms oi Sa (a = 1 . . .n) and S^^^ . One can use the set of equation (9) 
and (10) in order to find the linear relations between amplitudes ^^"^^ and 
According to Eq. (1) it will uniquely determine the matrix elements 
of Sa- The result reads: 

^ = f + £'5(o) L^f 

1 - f'5(o) 

Here r (r') and t (r) are block diagonal matrices with elements fmn,ai3 = 
5ali{ra)mn and tmn,Q/3 = 5ai3{ia)mn- The denominator of Eq. (11) describes 
the multiple reflection due to contacts within the dot. It can be verified, 
that 5* is unitary, S"^ S = 1, provided S^^^ and are also unitary. The 
analogous construction of ^-matrix can be in principle implemented for 
the more complicated layouts. 

To evaluate FCS one substitutes S into Eqs. (3) and (4) in the form 
(11). The matrices Sa are fixed by the choice of the properties of the con- 
tacts. However, the matrix S^^^ is random describing chaotic scattering 
inside the dot. This means that the answer should be averaged over all pos- 
sible realization of S^^^ . A reasonable assumption is that 5^°^ is uniformly 
distributed in the space of all unitary matrices [30]. 

Thereby one solves the problem for a single node. For less trivial system, 
comprising two and more nodes, the expression for the 5-matrix, corre- 
sponding to Eq. (11), becomes more involved. Moreover one has to average 
over unitary matrices separately for each node. Besides, to describe a sim- 
ple system of diffusive wire one has eventually to go to the limit of infinitely 
many nodes. All this makes the scattering approach extremely inconvenient. 
Fortunately, one can treat these problems with a more flexible semiclassical 
Green function method, which we outline in the next section. It is applica- 
ble, since we assume that the conductance of the system G Gq. 
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3. Circuit theory of FCS in multi-terminal circuits 

In this section we formulate the semiclassical theory of FCS of non- interacting 
electrons in case of multi-terminal geometry. Our approach is based on the 
Keldysh Green function method. We develop the scheme to evaluate the 
FCS in the arbitrary multi-terminal mesoscopic system. For that we will 
use the circuit theory of mesoscopic transport [19]. Next, we illustrate this 
scheme, considering the big fluctuations of current in the three-terminal 
chaotic quantum dot. In the end of the section we discuss the shot-noise 
correlations and give the convenient expressions that depend on the scatter- 
ing properties of connectors only, and do not involve the scattering inside 
the cavity. 

3.1. CIRCUIT THEORY APPROACH TO THE FCS. 

We start by introducing current operators /j, each being associated with 
the current to/from a certain terminal i. Extending the method of [18] we 
introduce a Keldysh-type Green function defined by 



Here we follow notations of a comprehensive review [25]: Xi ^.re time- 
dependent parameters, f3 is a 2 X 2 matrix in Keldysh space, H is the 
one-particle Hamiltonian that incorporates all information about the sys- 
tem layout, including boundaries, defects and all kinds of elastic scattering. 
We use "hat", "bar" and "check" to denote operators in coordinate space, 
matrices in Keldysh space and operators in direct product of these spaces 
respectively. The current operator in Eq. (12), associated with the terminal 
i, is defined as /j = J (i^x^t(p/m)^' VFj(a;). Here ^' is the usual Fermi 
field operator and VFj(x) is chosen such a way that the spatial integration 
is restricted to the cross-section and yields the total current through the 
given lead. The Eq. (12) defines the Green function unambiguously pro- 
vided boundary conditions are satisfied: G{t,t') = G{x,x';t,t') approaches 
the common equilibrium Keldysh Green functions Gf'\t — t') provided x,x' 

are sufficiently far in the terminal i. These G^^^ (t — t') incorporate informa- 
tion about the state of the terminals: their voltages Vi and temperatures 



In the following we will operate with the cumulant generating function 
("action") S{{xi}) defined as a sum of all closed diagrams 



(12) 



Ti. 




(13) 
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Here T-j- (T^) denotes the (anti) time ordering operator. One can see by tra- 
ditional diagrammatic methods [25] that the expansion of S{{xi}) in powers 
of Xii't) generates all possible irreducible diagrams for higher order corre- 
lators of Ii{t) and thereby incorporates all the information about statistics 
of charge transfer. If we limit our attention to low-frequency limit of cur- 
rent correlations, we can keep time-independent Xi- In this case, the Green 
functions are functions of time difference only and the Eq.(12) separates 
in energy representation. Then the action S{{xi}) can be conveniently ex- 
pressed via the average x-dependent currents Ii{{Xi\) in the following way 

i|| = /.({x<})-/|^(W.(5(.)) (14) 

where to denotes the time of measurement. Thus defined cumulant gen- 
erating functions allows to evaluate the probability for Ni electrons to be 
transferred to the terminal i during time interval to in accordance with the 
general relation (2). 

Up to the moment the above Eqs. (12) and (13) can be used to define the 
statistics of any quantum mechanical variable. However we are interested 
in the statistics of the charge transfer. Since the charge is the conserved 
quantity, the proper construction of current operators Ij requires the gauge 
invariance of the Hamiltonian. Therefore fully quantum mechanical scheme 
includes the "counting" fields Xi(t) as gauge fields. (See the contribution 
of L.Levitov in this book.) It means, that the initial physical Hamiltonian 
i?(q, p) of the system should be replaced by the x-dcpcndcnt Hamiltonian 
H-^ = iJ(q, p — ^Ts^jXjVFi). Thus the appearance of the "counting" 
fields in the problem is similar to the inclusion of a vector potential A(r) 
of an ordinary electromagnetic gauge field. The crucial difference is the 
change of sign of interaction at the forward and backward branches of the 
Keldysh contour, that is reflected by the presence of fs matrix. Then the 
X-dependent Green function obeys the equation of motion, written with 
the use of x-dependent Hamiltonian i?^. Accordingly, the "action" S{{xi}) 
is defined as a sum of all closed diagrams with respect to the interaction 
Hint = — H. 

Thus defined constructions are locally gauge invariant with respect to 
rotation in the Keldysh space. In particular, the average current (14) is 
a conserved quantity, provided the current operator is defined as /j = 
dH^{t)/dxi- Therefore it is possible to perform the local gauge transforma- 
tion ip' = exp{ — (i/2)f3 Xi-^i} V' in order to eliminate the x-dcpcndent 
terms from the equation of motion(12) [18, 20]. The x-depcndcncc of G is 
thereby transferred to the boundary conditions: the gauged Green function 
far in each terminal shall approach Gi{e) defined as 

Gi{e) = exp(ixir3/2)Gf (e) exp(-ixif3/2) (15) 
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Here (e) corresponds to the equilibrium Keldysh Green function suffi- 
ciently far in the given terminal i. The precise form of the Green functions 
Gf'\e) will be explicitly given below in the text. [See Eq. (19)] 

Let us also note that the modification of current operators !{ by the 
"counting" fields can be safely omitted at energies near the Fermi surface, 
where one can linearize the electron spectrum. This is possible in the semi- 
classical approximation when the variation of V-Fj is small at the scale of 
the Fermi wave length. Expanding the on the Fermi shell, one arrives to 
Eqs. (12) and (13). Afterwards it is feasible to verify the result (15), with 
making use of the semiclassical Eilenberger equations. (For the details we 
refer the reader to the paper of W.Belzig in this book.) 

In the present form, the Eq. (12) with relations (15,14) solves the prob- 
lem of determination of the FCS for any arbitrary system layout: one just 
has to find exact quantum-mechanical solution of a Green function prob- 
lem. This is hardly constructive, and we proceed further by deriving a 
simplified semiclassical approach. First, we note that even in its exact 
quantum-mechanical form the Eq.(12) possesses an important property. 
We consider the quantity defined similar to standard definition of current 
density, j°(x,e) = limj-^a,/ ( V'" — V^)G{x,x';e)/m. By virtue of Eq.(12) 
this quantity conserves so that 

df{x,e)/dx'' = (16) 

This relation looks like the conservation law of particle current at a given 
energy. However, this relation contains more information since it is a con- 
servation law for a 2 X 2 matrix current. 

Next, we construct a theory which makes use of this conservation law. 
We concentrate on the semiclassical Green function in coinciding points, 
G{x, e) = iG{x, x'; e)/Tru, where u is a density of states at Fermi energy. So 
defined Green function has been introduced in several semiclassical theories. 
[25, 26, 27] It satisfies the normalization condition G^ = 1. We relate the 
"current density" j to gradients and/or changes of G{x), very much like 
the electric current density is related to the voltage in circuit theory of 
electric conductance. Following the approach of the circuit theory [19], we 
separate a mesoscopic layout into elements: nodes and connectors, so that 
the G{x) is constant across the nodes and drops across the connectors. 
One may associate a graph with each circuit, so that its lines would 
denote the connectors, and internal and external vertices correspond to the 
nodes and terminals, respectively. (See Fig. 3) This separation of actual 
layout is rather heuristic, similar to separation of an electric conductor 
of a complicated geometry onto nodes and circuit theory elements. The 
bigger the number and the finer the mesh of the nodes and connectors, the 
better the circuit theory approximates the actual layout. The idea for this 
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Figure 3. The graph of the circuit theory, associated with a 3 terminal mesoscopic 
system. Gi, G2, Ga in the terminal are fixed by the boundary condition (15). T^-' define 
the transport properties of a connector. lik and Ijk. denote the currents, flowing from the 
nodes i and j into the node k. xi, X2, X3 are counting fields in the terminals. 



separation is completely analogous to the one considered in the section 2. 
The nodes are similar to the terminals, the difference is that G is fixed in 
the terminals and yet to be determined in the nodes. The G in nodes are 
determined from Kirchoff rules reflecting the conservation law (16): sum of 
the matrix currents from the node over all connectors should equal zero at 
each energy. For this, we should be able to express the matrix current via 
each connector as a function of two matrices Gij at its ends. 

To accomplish this task, consider the connector linking to nodes 

i and j, It can be quite generally characterized by a set of transmission 
eigenvalues Tn '^\l9, 27]. The problem to solve is to express matrix current 
lij via the connector in terms of GjQ) . This problem shall be addressed by 
a more microscopic approach and was solved in [27] for Keldysh-Nambu 
matrix structure of G. It is a good news that the derivation made in [27] 
does not depend on concrete matrix structure and can be used for the 
present problem without any modification yielding 



I = J- y / dE — [^"^^'J — . (17) 



Each connector (i,j) in the layout contributes to the total Xi-dependent 
action (14). The corresponding Sij contribution to the action should be 
found from the relation (14) and reads: [20] 



^^iW = ^E / ^^Trln [l + ^Tfe) ({G„G,}-2) 



(18) 



Now we are ready to present a set of circuit theory rules that enables 
us to evaluate the FCS for an arbitrary mesoscopic layout. 

(i) The layout is separated onto terminals, nodes, and connectors. 
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(ii) The Gj in each terminal j is fixed by relation (15) thus incorporating 
information about voltage, temperature and counting field % in each node. 

(iii) For each node k, the matrix current conservation yields a Kirchoff 

equation Zj^ = 0, where the summation is going over all connectors {i, k) 
attached to node k, and lik arc expressed with (17) in terms of Gj(fc). 

(iv) The solution of resulting equations with condition = 1 fixes Gk 
in each node. 

(v) The total action S{x) is obtained by summing up the contribu- 
tions Sij{{xi}) of individual connectors, those are given by (18): S{{xi}) = 
J2{i,j)Siji{Xi}) 

(vi) The statistics of electron transfer is obtained from the relation (2) . 
In the end of this subsection we discuss the limits of applicability of the 

whole scheme. By virtue of semiclassical approach, the mesoscopic fluctua- 
tions coming from interference of electrons penetrating different connectors 
are disregarded. So that, we assume that conductivities of all connectors 
are much bigger than conductance quantum e? /ttH. The same condition 
provides the absence of Coulomb blockade effects in the system. Besides, 
we have disregarded the possible processes of inelastic relaxation in the 
system. The latter can be eventually taken into account, since the use of 
Keldysh Green functions technique allows for perturbation treatment of 
interaction and relaxation. However, it would considerably complicate the 
scheme. The point is that the inelastic scattering would mix up the G{e) 
at different energies, so that one can not solve the circuit theory equations 
separately at each energy. 

3.2. THE STATISTICS OF CHARGE TRANSFER IN CHAOTIC QUANTUM 
DOTS 

As an illustration of the above scheme, we will consider in the second part 
of this section the FCS of the 3-terminal chaotic quantum dot. The system 
is sketched in the inset of Fig. 2. The heuristic circuit, associated with this 
mesoscopic system is shown by dashed lines. It includes only 3 terminals, 
3 arbitrary connectors, associated with the contacts, and the node {4}, 
representing the quantum dot itself. This separation is valid provided the 
cavity is in the quantum chaotic regime. (See [28] for definition). This regime 
corresponds to full isotropization of the Green function G{x,x',€) within 
the dot, so that G4{e) can be regarded as a constant at a given energy. 

Since the normalization G\ = 1 holds for each vertex, we use the 
parametrization Gk = gfc • t, gk ' Sk = 1- Here is a 3-D vector, and 
T = (Ti,f2,f3). With the use of this parametrization the anticommutator 
{Gi,Gk} is proportional to the unity matrix and takes the form of scalar 
product ^{Gi,Gk} = gj • gfc. In the absence of counting fields the Green 
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functions in the terminals corresponds to the equihbrium Keldysh Green 
function [25] 

" V -2(1 - fk) - 1 ; ' ^'''> 

where Fermi distribution function fk{E) = {exp[(£' — eVfc)/Tfc] + ac- 
counts for the bias voltages Vj- and the temperatures T/j in the terminals. 
The Xi-dependence of Gfc(x) is then given by Eq. (15). 

Wc see that Green function G/^ (x) = g4 • r in the dot is in fact the only 
function to find. For that, we proceed by applying the current conservation 
law, Yl\^i Ik,4 = 0, inside the dot. We present the currents I^^^ given by (17) 
in the form ^ 

4,4 = ■^Zk{gk ■ g4)[Gk,G4], (20) 

where the scalar function Zk{x) incorporates the information about trans- 
mission eigenvalues in each connector k: 

Zk{x) ^ ^i'''V[2 + T(t''\x - 1)]. (21) 

n 

It can be evaluated for any particular distribution p{T) of transmission 
eigenvalues in the given connector and completely defines its scattering 
properties. For a example, if we denote by Rq = irh/e^ the resistance 
quantum, then R^^ = 2RQ^Zf;{l) is an inverse resistance of the connector. 
One can also express the Fano factor = (T(l — T))/{T), associated with 
the given connector as F = 1 — 2{d/dx)log Zk{x)\^^^. 

With the use of Zk{x) the conservation law can be efficiently rewritten 

3 

as [X^ p^GkjGi] = 0, where = Z^ (gfc • g4). This relation suggests to 
k=l 

3 

look for the vector g4 in the form g4 = J2 p'^ Sk, with M(x) being 

k=l 

an unknown normalization constant. Using the normalization condition 
g4 • g4 = 1 we obtain the following set of equations 

3 3 

p' = Z,(^M-'Y,9ij{x)P^), M^=Y,gi^{x)py (22) 

i=i «,i=i 

where the scalar product gij{x) = Eiix) ' Sjix) between terminal Green 
function is expressed in terms of Fermi distributions as follows 

5.,(x) = (1 - 2/0(1 - 2/,) + 2e^fe->^^)/,(l - /,) + 2 e-^te->^^)/,(l - /,) 

The Green function G4 then is found from the solution {p*(x); (x)} of 
this set of equations. In the general situation the function Zk{x) takes the 
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form (21) and therefore Eq. (22) represents the set of non-linear equations. 
However, their solution can be relatively easy found numerically using the 
method of subsequent iterations. 

The total action can be found by applying the rule (v) of circuit theory 
and reads ^ 

S(x) = *ji2 J deS,{g,-g^) (23) 
1=1 

where 

3 

g, • g4 = M-\x) E 9ij{x)p\x) (24) 
j=i 

Here partial contributions Sk{x) from each connector in Eq. (23) has to be 
determined from the relation -^Sk{x) = —Zk(x), provided Sk{l) = 0. It 
follows from the Eqs. (14), (20) and normalization condition g4 • g4 = 1. 

We consider three particular types of connectors: tunncI(T), ballistic(B) 
and diffusive(D). Their corresponding contributions to action read as [18]: 

St{x) = -^(i?o/i?)(x-l), (25) 
Ssix) = -iRo/R)log[{l + x)/2], (26) 

Snix) = ~{R^/R)log\x+Vx^-l) (27) 

with R being a resistance of the connector. The tunnel connector represents 
the tunnel junction, so that T„ <C 1 for all n. Ballistic connector corresponds 
to the quantum point contact, with N open channels. The last expression 
comes from universal transmission distribution p{T) = Rq/2RT\/1 — T for 
any diffusive contact. 

Analytical results for FCS (23) can be readily obtained only for the 
system with tunnel connectors. To assess general situation we found g4 
for given Xi numerically. To find the probability distribution, we evaluated 
the integral (2) in the saddle point approximation, assuming Xi has to be 
complex numbers. Saddle point approximation is always valid in the low 
frequency limit we consider, since in this case both action S and number 
of transmitted particles iVj = lito/e ^ 1. Due to the current conservation 
law only two of three counting fields Xi ai's independent, and one can set 
X3 = 0. The relevant saddle point of the function ft{x) = 5'(x)+^XiAio/e+ 
iX2-^2io/c always corresponds to purely imaginary numbers {Xi)X2}- '^^^ 
probability reads P(/i, 12) ~ exp[— Q(x*)]. Evidently, ^{x*) is the Legendre 
transform of the action, and it can be regarded as implicit function on I{x*)- 

In the following we assume the shot noise regime eV ^ kT when the 
thermal fluctuations can be disregarded. The energy integration in (23) be- 
comes trivial, since /j(e) = or 1, and it is sufficient to consider only the 
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IjR/V 



Figure 4- The logarithm of the current probabihties in the 3-terminal chaotic quantum 
dot as a function of Is, under condition Ii — l2- The insert presents the system configura- 
tion. The resistances R of all connectors are assumed to be equal. 1 - tunnel connectors, 
2 - diffusive connectors, 3 - ballistic connectors. 



case Vi = V2 = 0, Vs = V . Any other possible setup can be reduced to the 
number of previous ones by appropriate subdividing a relevant energy strip. 
The results of these calculations are shown in Fig. 4 and 5. We see that the 
maximum of probability occurs at Ii = I2 = —V/3R, I3 = 2V/3R. This 
simply reflects the usual Kirchoff rules. The current distribution P(/i,/2) 
for a ballistic system is bounded. This stems from the fact, that each ballis- 
tic contact has a limited number of open channels. Therefore any plausi- 
ble current fluctuation in the given terminal may not exceed the threshold 
value of a current Ig = NV/Rq via the ballistic connector. Contrary to 
that, for the tunnel and diffusive type configurations any connector has 
an infinite number of partially open transmission channels. Therefore the 
range of current fluctuations is not bounded in this case. 

From Fig. 4 and 5 it follows that the relative probabilities of big current 
fluctuations increase in the sequence ballistic— >diffusive—>tunnel. To reveal 
the origin of this behavior we proceed by considering the shot-noise cross- 
correlations Pij, defined by Exp. (5). At zero frequency they can be found 

Taking into account. 



analytically via the relation Pij = — f - a a S 



«J to dxidx. 



that the action S{{xi}), regarded as function of Xi, is implicitly determined 
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Figure 5. The contour maps of the current distribution log[P(/i, 12)] in the 3-terminal 
chaotic quantum dot for different configurations of connectors, (a) - ballistic connectors, 
(b) - diffusive connectors, (c) - tunnel connectors. 



by Eq. (23) via the set of equations (22), we arrive to the following result 
P^j = -G,j{T + 64) + R-^aiaj{@ - 6^ - 6^) + dijR'^aiQi (28) 
Here = R/R,, R-^ = Z^ 

G^j = R-\aiaj - 6ija^) (29) 
is a conductance matrix of the dot, 

3 

94 = / dsf{l - /), Qk = FkJ deifk - /)^ Q = Y1 (30) 

1=1 



18 



DMITRI A. BAGRETS, YULI V. NAZAROV 



Ffc is a Fano factor of the k-th connector and /(e) = Yli=i (^ifii^) is the 
non-equhbrium distribution function within the dot. 

In case of balhstic contacts = and the first two terms of the re- 
sult (28) reproduce the expression for the noise power obtained with the 
use of "minimal correlation" approach[21]. At zero temperature it also co- 
incides with the result of random matrix theory [29] , as wc expected from 
the our consideration of scattering approach. At equilibrium Pij = —2TGij 
in accordance with the fluctuation-dissipation theorem. In the general situ- 
ation each non-ideal connector with Fj. ^ gives the additive contribution 
to the cross-correlation function, which is hnear in i^^.. Thus we conclude 
that the current-current correlations in the chaotic quantum dot result from 
two contributions. The first one, proportional to (T + 64), corresponds to 
the cross-correlation function P^j of the chaotic cavity with ideal ballistic 
point contacts, which stems from the fluctuation of the distribution func- 
tion / within the dot. The second contribution, proportional to Gj, reflects 
the noise due to connectors. 

To present the result for the cross-correlations it is also useful to in- 
troduce the (3 X 3) matrix F with elements Fij = Pij/el^.^ where Iy. = 
X]i=i '^^^ matrix F is a generalization of the Fano factor for the multi- 
terminal system. It is symmetric and obeys the relation Yl\^i Fik = 0, 
which follows from the current conservation law. At zero temperature for 
the symmetric setup, shown in Fig. 4, it reads as 

^ f 4 + 3F -2 -(2 + 3F) \ 

F = —\ -2 4 + 3F -(2 + 3F) (31) 

"^^ \ -(2 + 3F) -(2 + 3F) 4 + 6F / 

For a diffusive wire = 1/3 and for a tunnel junction F^ = 1. There- 
fore at fixed average currents through connectors the Gaussian's currents 
fluctuations will increase in the sequence ballistic— >-diffusive—>-tunnel. As it 
was mention before similar behavior is also traced in the regime of the big 
current fluctuations. The essential point here is that the cross-correlations 
always persist regardless the concrete construction of the connectors. 

4. FCS of Charge Transfer in Coulomb Blockade Systems 

In the preceding sections wc have considered the FCS of non-interacting 
electrons. It was assumed that the conductance G of the system is much 
greater than the quantum conductance Gq = e^/h. It is known, that un- 
der opposite condition, G < Gq, the effects of Coulomb blockade become 
important. This motivates us to study the statistics of charge transfer in 
the mesoscopic systems, placed in the limit of strong Coulomb interaction, 
G <C Gq. The electrons dynamics in this Coulomb blockade limit is fortu- 
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nately relatively simple, since the evolution of the system is governed by a 
master equation. The charge transfer is thus a classical stochastic process 
rather than the quantum mechanical one. Nevertheless the FCS is by no 
means trivial. In this section we elaborate the general approach to FCS in 
the systems, governed by master equation. 

We begin this section by presenting the general model, which dynamics 
obeys the master equation. Further on we proceed with the proof of the 
central result of this section for the FCS in the master equation approach. 
We illustrate the general scheme by considering the statistics of charge 
transfer through the Coulomb blockade island with 3 leads attached and 
compare the FCS in this case with the results of the preceding section, i.e. 
for non-interacting electrons. 

4.1. THE GENERAL MODEL 

The dynamics of various systems can be described by master equation. For 
our purposes it is convenient to write it in the matrix form: 

||p(t)) = -L|p(t)) (32) 

where each element Pnit) of the vector \p{t)) is the probability to find the 
system in the state n. The matrix elements of operator L are given by 

Here r„^„j stands for the transition rate from the state m to the state n, 
7„ presents the total transition rate from the state n. The L operator thus 
defined always has a zero eigenvalue, the corresponding eigenvector being 
the stationary solution of the master equation. 

Coulomb blockade mesoscopic systems always obey Eq. (32). The main 
advantage of the master equation approach is a possibility of non-perturba- 
tive treatment of the interaction effects. Below we first consider the master 
equation description of the general model system. This will prepare us to 
the next section where we derive the FCS method. 

The possible physical realization of the general model includes an ar- 
ray of Coulomb blockade quantum dots and a mesoscopic system with a 
number of resonant levels. Like in the preceding sections, it is convenient 
schematically to present the system as a graph (see Fig. 6) with each node 
a corresponding either to a single dot, a single resonant level or an external 
terminal. The line k = {a, (3), connecting the nodes a and f3, is associated 
with a possible electron transfer. Let M be the total number of nodes in this 
graph and L is a total number of lines. For a many-dot systems each line k 



(33) 
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1 n„ 




2 

Figure 6. The graph of the general model (See the main text). The terminals are 
connected with the system via external junctions 1,2 and 3. The nodes a, f3 and 7 are 
either resonant levels or dots, linked with each other by internal junctions fc's. The arrows 
denote the conventional direction of a current through each junction. 



corresponds to the tunnel junction. For systems with many resonant levels 
it corresponds to the possible transition between different levels, so that it 
does not necessary correspond to electron transfer in space. There are N 
external junction k = 1 . . . N, (N < L), connecting the terminals with the 
system. The currents through these junctions are directly measurable and 
hence are of our interest. 

The macro- or microscopic state of the general model is given by a set 
of occupation numbers \n) = \ni, . . . ^um)', na is equal to any integer for 
the array of quantum dots and refers to the excess charge on the island 
a; in case of many resonant level denotes the occupation number of 
a given level. Owing to the fact that Yin ^nm = 0, the L operator has a 
zero eigenvalue. There are the right, |po)) and the left, (go|, eigenvectors 
corresponding to this zero eigenvalue 

L\pq) = 0, (go|i = (34) 

We assume that they are unique. This means that the system does not get 
stuck in any metastable state. The vector \pq) gives the steady probability 
distribution and (go| = (1, 1, • • • , !)• 

It is also useful to present L operator in the form 

L = 7-f, f = X;(fi"^^+fi-)) (35) 

k=l 

where 7 is the diagonal operator in the basis |n) of the system configuration 

and f is associated with the electron transfers through the line k = 
{a J): 

7 = J^|n)7(n)(n|, f = ^ |n')r(±)(n)(n| (36) 
{n} {n] 
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The state \n') = |ni, . . . , n^, . . . , n^, . . . , um) results from the state \n) by 
appropriate changing the corresponding occupation numbers: = no, — (7^, 
n'p = nfi + where crfc = ±1 denotes the direction of the transition. 

4.2. THE FCS IN THE MASTER EQUATION 

In this section we derive the central result for the FCS of the charge transfer 
in the system, which dynamics obeys the master equation. We will solve 
this problem by making use of the property of the system, that its random 
evolution in time is the Markov stochastic process. 

In what follows we will partially use notations of the book [31]. Let us 
consider the time interval [— T/2,T/2]. Suppose the system undergoes s 
transitions at random time moments r^, so that 

+r/2 >Ti>T2> ■■■>Ts-l>Ts> -r/2 (37) 

This gives an elementary random sample Qg = {j\,k\,ai; . . .;Ts-,ks-,(Ts)- It 
corresponds to the set of subsequent events, when at time the tunneling 
happens via the junction fcj, (jj = ±1 being the direction of the transition. 
The samples C,s constitute the set of all possible random samples. 

Then one defines the measure (or the probability) dix{C,) at the set For 
this purpose we may very generally introduce the sequence of non-negative 
probabilities Qs{{Ti, h, ^i}) = Q{ti, ki,ai; . . . ;Ts, kg, cFs) > defined in Q 
so that 

-l-CXD 

dfi{C) = Qo + XI X ^«({'^«' ...dTs (38) 

s=l {ki,ai} 

The functions Q are normalized according to the condition 

/ dfi{0 = Qo+ 
Jn 

+00 r r ^ 

XX \ -\ QsiXriM.<yi\)^dT, = \ (39) 

T/2>Ti>->Ts>-T/2 

Each term in Exp. (38) corresponds to the probability of an elementary 
sample (,3- 

To accomplish the preliminaries, we remind the concept of a stochas- 
tic process. Mathematically speaking, it can be any integrable function 
A{t) = A{t, C) defined at the set CI and parametrically depending on time. 
It is sometimes convenient to omit the explicit ( dependence. We will use a 
" check" in this case to stress that the quantity in question is a random vari- 
able. Each stochastic process A{t, C) generates the sequence of time depen- 
dent functions < ^0(^)5 ^1(^5 ti, ki,ai), ... , As{t, {tj, ki, CTj}) [•. Its average 
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{A{t))Q over the space Q is defined as 

„ +00 
{Ait))^ = 1^ Ait, C)dM(C) = Aoit)Qo + J2 Yl (40) 

s=l {ki,(Ji} 

s 

As{t, {n, ki, ai})Qs{{Ti,ki, ai}) dn 

T/2>Ti>->Ts>-T/2 

The analogous expression should be used, for instance, to define the corre- 
lations {A(ti)B{t2))fi between any two stochastic processes. 

For the subsequent analysis we define the random process I^'^\t), cor- 
responding to the classical current through the external junction k < N: 

s 

i=l 

Here dm is included to take into account the direction of the jump and 
S{k — ki) = 5k,ki is the Kronecker 5 symbol. Given this definition at hand, 
we introduce the generating functional <S'[{xi(t)}] depending on N counting 
fields Xi(r), each of them associated with a given terminal i: 

eM-SiiXiit)}]) = ( exp{z ^ / drxn(r)/(")(r)/e} ) (42) 

with the average defined by Eq.(40). Let us note the remarkable similarity 

of this classical expression with the quantum mechanical action (13). As 
before, in the low-frequency limit of current correlations one may use the 
time-independent counting fields Xi- In this case the action S'dxi}] can be 
used to find the probability (2) of Ni electrons to be transferred through 
the terminal i during the time interval T. 

The above definitions were rather general than constructive, since the 
probabilities Q have not been specified so far. To proceed, one has to 
relate them to transition rates of the master equation. We assume that 
at initial time t = —T/2 the system was in the state {n^^^}. Then ran- 
dom sample (g determines the evolution of charge configuration {n^^^} 
l^(s-i)} . . . {n*^^'} — {n^^^} for subsequent moments of time. The choice of 
Cs specifies that the transition between neighboring charge states {riM^} and 
occurs at time Tj via the junction ki = (a,-,/?;). Therefore the se- 
quence {n*^*^} is given by the relation ni*. = rial ~ '^kn n^, = n^^j 
and for all 'J ^ ai and To determine the probability 
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Qsi{Ti, ki, (Tj}) we note that (i) the sample constitutes the Markov chain 
(ii) the conditional probability of the system to remain at state n^*^ between 
the times Ti^i and Tj is proportional to cxp[— 7(n(*))(Tj — Ti^i)]; (iii) the 
probability that the transition occurs via the junction ki during the time 
interval drj at the moment Tj is given by r^j^,^\n^^^)dTi. These arguments 
suggest that Q's have the form 

Qo = Vexp[-7(n('^)r] (43) 
Qs{{Ti,h,a,}) = Z,' cM-lin^''^)iT/2 - n)]Ti['\nW) 

exp[-7(n(^))(ri - r2)]rg^^(n(2)) . . . exp[-7(n(^-i)) 

(r._i - r.)]r(-)(nW) exp[-7(nW)(r, + T/2)] 

where the constant Zq should be found from the normalization condi- 
tion (39). As we will see below, Zq = 1. 

The above correspondence between the random Markov chain and the 
probabilities Q's (43) allows one to evaluate the generating function (42). 
By definition (41) for any given we have 




exp<i^ / dTXn{r)I^'^\r,(s)/e> ='[{exp{iaiXki{Ti)} 



It is assumed here that Xki = if the transition occurs via internal junction, 
ki > N, thus no physically measurable current is generated in this case. 
The averaging of the latter expression over all possible configurations CI 
with the weight dn{C) yields 

Z[{x^{T)}] = exp(-5[{x.(r)}]) = ^0+ (44) 

+00 „ „ s 

EE y Qi{{n,ki,a,})l[dn 

s=l{ki,ai} T/2>Ti>->Ts>-T/2 

The resulting expression resembles the normalization condition (39). Here 
the x-dependent functions Q^{{Ti,ki,(Ti}) arc defined similar to probabil- 
ities (43) with the only crucial difference that the rates T^i^\n) should be 

replaced by r['^'*(n) exp{icjfc Xfc(''"fc)} if A; < iV. 

The expression (44) can be written in the more compact and elegant 
way. For that, we introduce the x-dependent linear operator defined as 

4(r) = 7-r;,(r), (45) 

N L 



TxW = V(fWe^X'=(-)+f(-)e->^'=«)+ V (f(+)+f(- 



A^\-k ^ ' -k ^ ; ' ^ ■■ - - ■ 

fc=l k=N+l 
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In line with consideration above we multiplied each operator Fj^ {k = 
1. . .N), that corresponds to the transition through the external junction, 
by an extra x-dependent factor e*^*^(^^. The diagonal part and internal tran- 
sition operators f with k > N remained unchanged. Then we consider 
the evolution operator U^{ti,t2) associated with (45). Since L^{t) is in 
general time-dependent, U^(ti,t2) is given by the time-ordered exponent 

U^{h,t2) = Trew[- J\j{t) - t^{r))dT] (46) 

The similar construction is widely used in quantum statistics. The differ- 
ence in the present case is that the operator f7^(ti,t2) at X = gives the 
evolution of probability rather than the amplitude of probability. 

With the use of evolution operator (46) the generating function (44) 
can be cast into the form 

Z[{x^{r)}] = {qo\UxiT/2, -T/2)|n,) (47) 

To prove it wc argue as follows. We exploite the fact that 7(t) and T{t) 
commute under the sign of time-ordering in Eq. (46) and regard r(r) as a 
perturbation. This gives the matrix element {qo\U~^(T/2,—T/2)\ns) in the 
form of series 

{qoPxiTf^, -T/2)\ns) = (go|e-^^|po>+ (48) 

+°° fT/2 

T/2>Ti>->Ts>-T/2 

s 

rl?^(n)e'''i^'=i(^'=i) . . . f^l'\Ts)e'''^''^^^^'^^\vo) n 

i=l 

It follows from the definition (43) that each term in this series corresponds 
to the function Qs ({t^, /cj, CTj}), namely 

Qo = (go|e"^^|ns) 

Q?({ri, Gi}) = {qo\Tr exp{- J^^^^^ 7(T)dr} (49) 

ri? Vi ) e'''' ^'=1 (^'=1 ) . . . f ) e'"'^"' Vs ) 

Therefore Exp. (48) and (47) are reduced to the previous result (44). 
This completes the proof. Note, that owing to the property (34), Zq = 




Multiterminal counting statistics 



25 



{qo\exjp(—TL)\ns) = 1 identically at x = 0. Therefore the probabilities 
(43) are correctly normalized. 

The Exp. (47) for the generating function Z[{xj(t)}] depends on the 
initial state \ns) of the system. It can be shown that the choice of \ns) 
does not affect the final results. We assume that Xk{t) — when t 
—T/2. Physically, it means that the measurement is limited in time. To be 
specific one may assume that Xki^) = when —T/2 < t < —T/2 + At and 
Xki't) 7^ if i > —T/2 + At. If the time interval At is sufficiently large as 
compared with the typical transition time T"^, then the system will reach 
the steady state during this period of time. The latter follows from the fact 
that exp {-LAt)\ns) \po) when At > T^^ Thus one can substitute In,) 
to \po) in Exp. (47). Assuming also the limit T oo, we arrive to the main 
result of this section 

+ 00 

exp(-5[{xi(t)}]) = (go|r,exp{- J 4(r)dr}|po) (50) 

—oo 

We see that the generating function can be written in the form of the 
averaged evolution operator. This operator corresponds to master equation 
with the rates modified by the counting fields Xi(T). 

Further simplification is valid in the low frequency limit of the current 
correlations, lo <^ T. (Here F is a typical transition rate in the system.) 
This allows to set Xk{t) = Xk when < t < to and Xkit) = otherwise, 
where to is a time of measurement. The action (50) then reduces to the 

S{{Xi}) = toAmin({xO) (51) 

where Ajnin({Xj}) is the minimum eigenvalue of the operator L^. Thus 
the problem of statistics in the Coulomb blockade regime, provided the 
transition rates in the system are known, is merely a problem of the linear 
algebra. 

4.3. THE STATISTICS OF CHARGE TRANSFER IN COULOMB 
BLOCKADE ISLAND. 

In this subsection we use the developed method to study the current statis- 
tics in the thrcc-tcrminal Coulomb blockade island. Its equivalent circuit 
is shown in Fig. 7. This circuit is an extension of the usual set-up of a 
conventional single electron transistor [32]. 

The island is in the Coulomb blockade regime, Rk 3> Rq = 2Trh/e^. In 
order to observe the Coulomb blockade effect the condition fc^T <^ Ec = 
e^/2Cs should be also satisfied. Here Ec is the charging energy of the island 
and Cs = Yl^=i Ck + Cg. We assume the temperature to be rather high. 
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Figure 7. The equivalent circuit of the three-terminal Coulomb blockade island. The 
voltages f/i{2) are used to control the bias between the 3d and the 1st (the 2nd) terminals: 
Ui^2) = V3 — Vi(2). The third terminal is biased at voltage V3 — {Ui + U2)/3 with respect 
to the ground. This setup assures the condition V1 + V2 + V3 = 0. Then gate voltage, Vg, 
can be used to control the offset charge qo = CgVg on the island. 



ksT S> AE, with AE being the mean level spacing in the island, so that 
the discreteness of the energy spectrum in the island is not important. We 
also disregard the possible effects of co-tunneling. 

Under the above conditions the 3-terminal island is described by the 
"orthodox" Coulomb blockade theory. In this theory the macroscopic state 
of Coulomb blockade island is uniquely determined by the excess charge 
Q = ne, which is quantized in terms of electron charge (— e). The charge 
Q can be changed only by ibe in course of one tunneling event. There- 
fore the master equation connects the given macroscopic state n with the 
neighboring states tt, it 1 only. The corresponding rates T^-i-]^^ — ^ of these 
transitions are equal to the sum of tunneling rates across all junctions: 

= ^k=i ^nii^n- ^^^6 tunnehng rate r^^^^„ across the junction k 

(k) 

can be expressed via the electrostatic energy difference AE^^J.-^^^^ between 
the initial (n) and final (nil) configurations 

e -ftfc 1 - eM-AE'^l-^^^/kBT] 

(k) 

The evaluation of AE^^J.-^^^,^ can be done along the same lines as in the case 
of single electron transistor [32]. 

According to the general result (51) of the preceding subsection, one can 
find the FCS of the charge transfer through the island, by evaluating the 
minimum eigenvalue Amin of the matrix L^. In case under consideration 
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this problem is reduced to the eigenvalue problem of the three-diagonal 
matrix: 

(A - 7„K + r^^„+lPn+l + r^^n-lPn-l = (53) 

where 7n = r„^n-i + T^^n+i, and r^^„±i = Ef=i rSi^„e±^>^'=. The 
index +(— ) corresponds to electron transition from (to) the island. 

To assess the FCS we have treated the related linear problem (53) nu- 
merically. We restrict the consideration to sufficiently low temperatures 
ksT ^ Ec, so that the temperature dependence in rates (52) is non- 
essential. In this case T'j^\n) = AE^J^l^^J {e^ Rk) when A£;2i^„ > 

and T^^\n) = otherwise. Wc have also assumed that U2 > Ui (See 
Fig. 7). The corresponding ^-dependent rates can be found from Exp. (52) 
and (45) and read as follows 

r^+i^„ = 4+)(n)e'>^3 + r(+)(n) e{q - 1/2 - n)e'^^ 
r^_i^„ = r5"V)e~'^' + ^i~\n) e(n-q- l/2)e-'^^ (54) 



where 



and 



rf)(„) = „«^(„ + ^^±lj_ (55) 



(+) _ CiUi + C2U2 (-) _ {Ci + C^)U2-CiUi 



a 



(±) _ j_ C2U2 - {C3 + C2)Ul 

The effective capacitances Ck are defined as Ck = Ck + Cg/3 and the point 
q is given by the relation 

eq{Ui , U2 , Vg) = C2U2- {C3 + C2)Ui- CgVg (56) 

The value q is non-integer in general. It satisfies the condition r^"-* {q + 
1/2) = T^^\q — 1/2) = 0. The dimension of the L^-matrix is equal to 
'T-max— 'T-mirn where nniax(™min) Can be found from the conditions Fg (n) > 
(r^^^(n) > 0). The value enmax, (e^^min) gives the maximum (minimum) 
charge that can be in the island for a given voltages Ui, U2 and Vg. 

We can see from the Exp. (54) that there are four elementary processes 
of charge transfer in the system at low temperatures, each of them being 
associated with the pre-factor e^*^*" . The factors e*^^ and e"*^^ correspond 
to the charge transfer from the third terminal into the island and from the 



28 



DMITRI A. BAGRETS, YULI V. NAZAROV 




Figure 8. The matrix F of auto- and cross- shot noise correlation versus voltage Ui for 
the 3-terminal quantum dot setup. Parameters are shown on the plot. (1) - Fn, (2) - 

\F^2l (3) - |Fi3|, (4) - F22, (5) - IF23I, (6) - ^^33. 




0.0 0.5 
Offset charge (q /c) 



1.0 



Figure 9. The matrix F of auto- and cross- shot noise correlation versus the offset 
charge for the 3-terminal quantum dot setup. Parameters are the same as on the Fig. 8. 
The voltage C/i = (72/4 = 1.25e/CE. (1) - -Fii, (2) - |fi2|, (3) - \Fri\, (4) - F22, (5) - |f23|, 
(6) - F33. 



island into the second terminal, respectively. Hence, the random current 
through the 3d (2nd) junctions always has the positive (negative) sign. 
Two factors e^^'^ stem from the charge transfer through the first junction 
in the direction either from the island into the first contact or vice versa. 
Therefore the current Ii fluctuates in both directions. 
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Figure 10. The logarithm of current distribution lnP(7i, I2) in the 3-terminal quantum 
dot as a function of current Ji, under condition I2 ~ {12)- Parameters are shown on the 
plot. (1) - Ui = 1.25 b/Ce, (2) - C/i = 2.0e/CE, (3) - = 4.0c/Ce, (4) - (7i = 10.0 b/Ce; 
curve (5) corresponds to the non-interacting regime. 



Let us consider the shot noise correlations in the system. In Fig. 8 we 
give the illustrative example of the voltage dependence of the shot noise 
correlations for a certain choice of parameters. The definition of matrix 
-Ffcm is the same as we have used in the end of section 3. The Coulomb 
blockade features are strongly pronounced for an asymmetric setup only. 
The results shown in Fig. 8 correspond to Ri = R3 = R2/W, Ci = C2 = 
C3 and U2/U1 = 4. In Fig. 9 we show the dependence of the shot noise 
correlations on the offset charge for the same set of parameters and the 
value of Ui = 1.25 e/Cs. The special points of both these dependences 
occur when either nmim ^^max or the integer part of q are changed by ±1. 
As the result one observes multi-periodic Coulomb blockade oscillations in 
the offset charge dependences. 

We now proceed with the evaluation of the FCS. The action S{{xi}) 
has been calculated with the use of (51). To find the probability distri- 
bution P we have evaluated the Fourier transform (2) in the saddle point 
approximation. It is applicable here, since we consider the low frequency 
limit only, a; <C F. In this limit both action S and number of transmit- 
ted particles Ni = lit^/e » 1. Due to the current conservation Ylk^k 
= 0, only two currents are independent and the action S{{xi\) depends 
on the differences Xij = Xi ~ Xj only. In what follows we have chosen Ii 
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(b) 

Figure 11. The contour maps of the current distribution log[P(/i, 72)] in the 3-terminal 
quantum dot. U\ = (72/4.0 = 1.25 c/Ce- (a) - Coulomb blockade island. Parameters are 
the same as in Fig. 10. (b) - Chaotic quantum dot. 



and I2 as the independent variables to plot the logarithm of probability 
lnP(/i,l2). In the saddle point approximation, with the exponential accu- 
racy, it is given by P(/i,l2) ~ e~^^^ \ Here x* is ^ saddle point of the 
function ^{x) = S{x) + ixihto/e + 2x2-^2^0/6. The results for lnP(/i,/2) 
are shown in Figs. 10 and 11(a). From the contour map on Fig. 11(a) we 
see that P(/i, I2) is non-zero in the region Ii < 0, I2 < and in the region 
-^1 < I -^2 1 provided /i > and /2 < 0. This range of plausible current fluctu- 
ations stems from the ^-dependence of rates (54). Any current fluctuation 
automatically satisfies the constrain Ylk-^k = and conditions /2 < and 

Before discussing the results, let us set the reference point for such dis- 
cussion. This reference will be the results of the previous section. We con- 
sider the FCS in the three-terminal chaotic quantum dot when its contacts 
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are tunnel junctions with resistances R^^ 3> e^/nh. In this hmit the effects 
of interaction are neghgible and electrons are scattered independently at 
different energies. Provided U2 > Ui, the generating function S{{xi}) in 
the given case is a sum of the two independent processes (23) 

'S'(Xl>X2,X3) = S'l(xi,X2,X3) + 'S'2(X1,X2,X3) (57) 

Here 

5i(xi,X2,X3) = ^{Ci + Gs + Ga- 
sj (Gi +G2- G^Y + 4G3e'X3(Gie-^xi + G^e-i^-^)^ 

^2(xi,X2,X3) = ^^^~f^'° {G, + G2 + G3- 
^ (Gi + Gs - G2Y + 4G2e-'X2(Gie*xi + G3e'X3)| 

and Gfe = R^^ are the conductances of the junctions. 

The logarithm of probability In Po(-^i) -^2)) evaluated with the use of 
statistics (57), is shown by the dashed line in Fig. 10. Its contour map for 
the same values of parameters is separately presented in Fig. 11(b). The 
maximum of lnPo(A)-^2)) as expected, occur at Ii = I2 = Ui/'iRi. 

Comparing the FCS in the Coulomb blockade and non-interacting limits 
we can draw the following conclusions. In spite of the different regimes, we 
see that the qualitative dependence of probabilities versus the currents is 
similar for both cases. The probability distribution in both cases has a 
single maximum, corresponding to the average values of currents. The tails 
of distribution are essentially non-Gaussian both in the weak and strong 
interacting limit. The statistics approaches to the Gaussian- type one in 
the strong Coulomb blockade limit only, when the applied voltage to the 
system is only few above the Coulomb blockade threshold. (See curves 1 
and 2 in Fig. 10) At higher applied voltages the probability distribution has 
a tendency to approach to the current distribution of the non-interacting 
system. However, they never become identical, even in the limit Ui^2 S> 
e/Gs. (Curves 4 and 5). The same is true for the shot-noise correlations. 
Generally, we conclude that the Coulomb interaction always suppresses the 
relative probabilities of big current fluctuations. This stems from the fact 
that any big current fluctuation in Coulomb blockade dot is related with 
the large accumulation (or depletion) of the charge on the island. The latter 
process results in the excess of electrostatic energy. Therefore the relative 
probability of such fluctuation is decreased, as compared to the probability 
of the similar current fluctuation in the non-interacting regime. 
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Figure 12. The single resonant level system, formed by the two tunnel barriers. The 
resonant level in the quantum well is shown by the dashed line. 



5. The equivalence of scattering and master equation approaches 
to the FCS 

In this section we evaluate the FCS in the generic case of a single resonant 
level, shown in Fig. 12. We consider only the non-interacting spinless elec- 
trons and demonstrate the equivalence of scattering and master equation 
approaches to the FCS in the framework of this model. 

We start by considering the FCS in the master equation framework. 
It is applicable, provided the applied voltage or the temperature are not 
too low, i.e. max{ey, /c^r} ^ hV^^^. Here ^l{R) denote the quantum- 
mechanical tunneling rates from the left (right) electrode onto the resonant 
level. The system can be found in the two microscopic states only: one 
with no electron on the level, and another with a single electron. Then the 
transition rates involved reads as 

Fi^o = TLfL{ei)+TRfR{e,) (58) 
Fo^i = Tl[1- fL{ei)\+TR[l- fR[ei)] 

Here the indices {0} and {1} denote the microscopic state with no and one 
electron on the level. Fermi function fL(R){^) = + fJ'L{R))/kT])~^ 
accounts for the filling factor in the left (right) lead and is the position 
of the resonant level. 

Following the definition (45) and the expression for the rates (58), the 
L^-matrix of the single resonant level model reads as 



(59) 
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where 

ri^o(x) = Ti/ie-^^i + TRfne-'^^ (60) 

ro^i(x) = ri(i - /L)e*>^i + r,^(i - /«)e*>^^ 

Evaluating the minimum eigenvalue of this matrix we obtain the current 
statistics in the following form 

5(x) = |{r^ + r«-V^} (61) 

VU) = (Xl + Tr^ + ATlTr [/(-)(e.)(e-^>^ - 1) + /(+)(e.)(e^^ " 1)] 

Here /(_)(€,) = /L(ei)[l - fR{ei)], /(+)(ei) = /R(ei)[l - fUei)] and x = 
Xi - X2- 

Since the electrons above are assumed to be non-interacting, one might 
have come to the same result in the framework of the pioneering approach 
by Levitov et.al.[12]. We will show now that it is indeed the case. 

According to Ref. [12] the general expression for the current statistics 
through a single contact reads as 

^(^) = -^E / ^eln{l + r„(6)x (62) 

n 

(/L(e)[l - fRie)]ie~'^ - 1) + /«(e)[l - hieW^ - 1)) } 

It is valid for any two-terminal geometry provided the region between two 
electrodes can be described by the one-particle scattering approach. r„(e) 
is a set of transmission eigenvalues which are in general energy-dependent. 
For a resonant level there is a single resonant transmission eigen- value Tr{e), 
its energy dependence being given by the Breit-Wigner formula 

The result (62, 63) is more general, than Exp. (61), obtained by means 
of master equation. If electrons do not interact, the Exp. (62) is valid for 
any temperature. It can be simplified in the regime ksT ^ HT. As we will 
show, in this limit the general result (62) coincides with the result (61) of 
master equation. It is easier to perform the calculation if one first evaluates 
the x-dependent current J(x) = {ie/to)dS/dx- It reads 

HX) = ^jde [/(+)(e)e*^ - /(-)(e)e-'^] x (64) 
[T~\e) + [/(_)(e)(e->^ - 1) + /(+)(6)(e^>^ - 1)] } 
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Let us consider the situation when the resonant level is placed between 
the chemical potentials ij,l{r} in the leads. Since we assumed that ksT » 
^L(R)7 the main contribution comes from the Lorcntz peak and one can put 
e = ej in the Fermi functions. Therefore we left only with the two poles 
ei(2) = ± iy^P(x)/2 under the integrand (64). Closing the integration 
contour in the upper or lower half-plane we arrive at 

/(X) = er^r^ [/(+)(e,)e^>^ - /(_)(e,)e->^] 

Integrating it over x one finds for the S{x) = (to/ie) I{x')dx' the re- 
sult (61) obtained by means of master equation. 

Thus, we have verified the correspondence between two approaches to 
statistics in the non-interacting regime. We have shown that one can re- 
produce the statistics (61) on substituting Tr{e) into the Exp. (62) and 
assuming the regime ksT ^ hT. As it was discussed previously, this is 
the condition, when the master equation approach, and hence its conse- 
quence (61), are valid. 

6. Summary 

We have reviewed here a constructive theory of counting statistics for elec- 
tron transfer in multi-terminal mesoscopic systems. We have covered two 
opposite limit of weak and strong interaction. For the case of weakly inter- 
acting electrons, when the conductance of the system G » Gq , the theory 
of FCS reduces to a circuit theory of 2 x 2 matrices associated with Keldysh 
Green functions. In the Coulomb blockade limit, G <^ Gq, the FCS meth- 
ods turns out to be an extension of the usual master equation approach. We 
have applied these methods to study the FCS of charge transfer through 
the three-terminal quantum dot. Surprisingly, the FCS has a similar qual- 
itative features both in weakly and strongly interacting regimes. We found 
that Coulomb interaction suppresses the relative probabilities of big current 
fluctuations in the dot. We have also reviewed the scattering approach to 
FCS in multi-terminal circuits. Then by considering the generic model of 
a single resonance level, we have established the equivalence of scattering 
and master equation approaches to FCS. 

The theories presented enables one for easy theoretical prediction of 
the FCS for a given practical layout. Thereby they facilitate experimental 
activities in this direction. Up to now, only the noise has been measured. 
In our opinion, the measurements of FCS can be easily performed with 
threshold detectors that produce a signal provided the current in a certain 
terminal exceeds the threshold value. If the threshold value exceeds the 
average current, the detector will be switched by this relatively improbable 
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fluctuation of the current. The signal rate wiU be thus proportional to the 
probability of these fluctuations P{I), the value given by the theory of FCS. 
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